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INTERACTION OF THE SPACE SHUTTLE 


CONTROL SYSTEM WITH POGO 

By A. Stewart Hopkins and William F. Davis 
McDonnell Douglas Astronautics Company 


SUMMARY 


Pogo and control instabilities have been of major concern in the design 
of aerospace vehicles; however, the symmetry of past vehicles and the 
associated absence of axial-lateral coupling have precluded the interaction 
of the pogo and control systems. 

Currently proposed Space Shuttle Systems consist of an orbiter riding in 
parallel with various combinations of propellant tanks and boosters. These 
asymmetric configurations result in coupled axial and pitch modes of vibra- 
tion. The pogo system, sensitive to axial motion, and the control system, 
sensitive to pitch motion, will interact because of the coupled modes. To 
investigate this phenomenon, detailed models of a representative Space Shut- 
tle structure and the associated feedline, control, and engine systems were 
acquired or developed. A model generated by coupling the subsystems in an 
appropriate manner was evaluated for stability by calculation of the complex 
eigenvalues. Eigenvalues of the coupled system were compared with those of 
systems with the engine or control system inactive. It was determined that 
even though a vehicle is stable with respect to both the pogo and control sys- 
tems individually, the coupled system may have a strong instability. To allow 
the traditional independent analysis and design of the pogo and control sys- 
tems, a method was developed for determining the stability of the coupled 
system in terms of the stability of the pogo and control systems separately. 
Although the method becomes imprecise for the case of very close roots, a 
simple extension provides upper and lower bounds for the stability. 



INTRODUCTION 


Oscillations due to coupled propulsion-structure instabilities, pogo, and 
control- structure instabilities have been experienced on numerous aerospace 
vehicles. The source of the pogo instability may be envisioned as follows. 
Propellant oscillations induced by axial structural oscillations produce varia- 
tions in engine inlet pressure and propellant flow. These result in thrust per- 
turbations which reinforce the structural motion. Similarly, control signal 

oscillations , induc'echby" late ral s true tural-osc illation s r .p rod u c e_ control, for ces 

which reinforce the structural motion. The axial and lateral modes are almost 
entirely uncoupled for traditional axisymmetric vehicles. Methods have been 
developed to evaluate pogo stability with respect to axial modes and control 
stability with respect to lateral modes. Requirements will exist to ensure the 
stability of the Space Shuttle in this conventional sense. 

The Space Shuttle configuration, however, introduces the possibility of a 
new type of instability produced by pogo-control interaction. Currently pro- 
posed Space Shuttle Systems, such as illustrated in fig. 1, exhibit a high 
degree of asymmetry, which results from coupling an orbiter with various 
combinations of boosters and propellant tanks. As has been confirmed by 
preliminary analysis, the structural asymmetry results in strongly coupled 
axial and lateral modes. Structural oscillations, therefore, induce responses 
in both the propulsion and control systems. Additionally, the long lateral 
feedlines characteristic of many Space Shuttle designs allow propellant oscil- 
lations to be induced by lateral structural motion. The combined regenerative 
forces may be of substantially greater magnitude and different phase, intro- 
ducing the possibility of instabilities that would not be predicted by conventional 
methods. 

Although the study conclusions are relevant to any Space Shuttle configu- 
ration exhibiting structural asymmetry with attendant axial-lateral coupling or 
lateral feedline runs, the McDonnell Douglas Corporation (MDC) "recom- 
mended Phase B fully reusable baseline configuration" was selected for 
detailed analysis. This representative configuration was selected so the 
existing structural idealization, control and engine system models, and feed- 
line design could be utilized. Analysis was limited to the pitch plane since 
symmetry precludes coupling with the remaining degrees of freedom, and to 
the LOX line because the fuel feedline runs are relatively short and the engine 
is insensitive to fuel side perturbations. The coupled system model was 
developed to accurately reflect the system behavior from 0. 07 to 7 Hz. This 
includes the dominant frequency ranges of the control system(0. 2 to 3 Hz), 
engine system (0. 1 Hz and up), and the fundamental feedline and structural 
modes (0.7 to 7 Hz). 

The purpose of this study was to determine the nature of the pogo-control 
system interaction and to develop evaluative stability analysis techniques. 

The stability of the coupled system was evaluated with closed- loop analysis 
procedures. These methods were found to be substantially more effective 
than open-loop or time-domain procedures. After the stability of a baseline 
configuration was evaluated at selected burn times, the period of maximum 
control participation, 25 percent burn, was selected for evaluating sensitivity 
to parametric variations of the structural, feedline, control, and engine 
systems . 


2 




SYMBOLS 


_-A 


a 



A 

FT 

A G 


A 


R 


B 


b 

C 


D 


d 

D < 

D 

Dc 

k. 

D 

E 

F 


EE 


SS 


control conditioned pitch rate 

real part of system eigenvalue 

area of feedline at engine 

area of feedline segment 

area of feedline at tank 
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effective area ratio defined in equation (4) 

control conditioned attitude error 

imaginary part of system eigenvalue 

control command signal 

aerodynamic lift coefficient gradient 

coupled system damping matrix 

feedline segment diameter 

control system damping matrix 

engine system damping matrix 

structure -feedline dampipg matrix 
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axial force at gimbal 
lateral force at gimbal 
moment at gimbal 
acceleration of gravity 
feedline modal force vector 

height of liquid surface above center of mass 
structure modal force vector 
identity matrix 

\TT 

engine moment of inertia x 12 
interaction coefficient 
coupled system stiffness matrix 
aerodynamic spring defined in equation (2) 
control system stiffness matrix 
engine system stiffness matrix 

a, P = R (for reaction), C (for cavitation), RC (for both 
reaction and cavitation), and F (for free); feedline stiff- 
ness matrix partitions defined by equations (10) and (11) 

structure-feedline stiffness matrix 

augmented structure stiffness matrix partition 

ullage spring rate 

control attitude gain 

control attitude rate gain 

feedline segment length 

center of pressure to center of mass distance 
coupled system mass matrix 
actuator moment 
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control system mass matrix 
engine mass x 12 
engine system mass matrix 
feedline mass matrix partition 
tank fluid mass 

a, P = S (for structure), F (for fluid); structure- 
tank mass matrix partitions defined by equation (6) 

mass flow from tank 

structure -feedline mass matrix 

slosh mass 

a, p = S (for structure), T (for tank fluid); augmented struc- 
ture mass matrix partitions defined by equations (7) and (8) 

Mach number 

ullage gas pressure 

suction pressure at engine inlet 

pressure at tank sump 

dynamic pressure 

engine fluid inertial force 

engine fluid inertial displacement 

feedline force vector 

feedline displacement vector 

feedline force at tank 

feedline displacement at tank 

total displacement tank fluid center of mass 

feedline reaction force vector 

feedline reaction displacement vector 

vector whose partitions are and Q^_ defined by 

equation (11) 
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vector whose partitions are q R and q„ defined by 
e quation (11) ^ 

vector whose partitions are q R and X defined by 
equation (12) 

augmented structure force vector 

augmented structure displacement vector 

slosh generalized displacement 

Laplace operator 

aerodynamic reference area 

engine thrust x 12 

feedline wall thickness 

tank ullage gas volume 

tank fluid volume 

weight flow through engine 

global state vector 

engine axial displacement 

displacement structural tank mode 

tank fluid center of mass motion due to flow 

axial displacement of tank bottom 

lateral displacement at engine center of mass 

lateral displacement at gimbal 

aerodynamic angle of attack 

LOX bulk modulus 

control gain parameter 

engine gimbal angle 

control frequency parameter 

feedline modal amplitude vector 

structure modal amplitude vector 
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ullage gas specific heat ratio 

gimbal to engine center of mass distance 

metric scale factor 

system natural frequency 

aerodynamic-frequency — 

feedline frequency diagonal matrix 
structural frequency diagonal matrix 
slosh frequency 
feedline mode shape matrix 

vector of structure rotational mode shapes at attitude 
reference unit 

structure mode shape matrix at reaction coordinates 

structure mode shape matrix at reaction coordinates 
and engine 

vector of structure rotational mode shapes at rate gyro 
vector of structure axial mode shapes at gimbal 

vector of structure lateral mode shapes at gimbal 

vector of structure rotational mode shapes at gimbal 

vector of structure mode shapes for fluid tank mode 

vector of structure mode shapes at tank bottom 
lateral mode shape in rigid -body plunge 

vector of structure lateral mode shapes at tank bottom 

rotational mode shape in rigid-body pitch 
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LOX mass density- 

pitch rotation at attitude reference unit 

engine rotation 

rotation at gimbal 

pitch rotation at rate gyro 

system fraction of critical damping 

aerodynamic fraction of critical damping 

system fraction of critical damping, control only 

feedline fraction of critical damping diagonal matrix 

system fraction of critical damping, pogo only 

system fraction of critical damping, pogo and control 

structure fraction of critical damping diagonal matrix 

slosh fraction of critical damping 

system fraction of critical damping, neither pogo nor 
control active 


ANALYTICAL MODELS 


To investigate the interaction of the pogo and control systems, the repre- 
sentative Phase B fully reusable baseline configuration was modeled in detail. 
The total system model was a composite of structural, fluid, control, and 
engine system analytical models. 

The available structural model (representing a somewhat shorter version 
of the Shuttle) was scaled to the baseline configuration length, aerodynamic 
and slosh effects were incorporated, and fluid motion into the tank was added 
to facilitate coupling to the top of the feedline. The feedline model was devel- 
oped from the feedline design and includes structural coordinates at bends for 
coupling to the structure, and weight flow for coupling to the engine model. 
Springs were incorporated to reflect pump cavitation and the design accumu- 
lator. The Phase B baseline control system was incorporated, including 
inertial effects of the gimbaled engine. A preliminary analytical engine model 
was obtained from an engine manufacturer. These four analytical models 
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were incorporated in a composite model by transforming each to a common 
set of variables and applying the appropriate compatibility and equilibrium 
relationships. The common set of variables includes pressures and flow as 
well as displacements and is therefore called a state vector. 


Structure 


The structural model was generated by making three enhancements to the 
available Phase B baseline idealization. Certain aerodynamic effects were 
included because of their potential effect on the control system. A slosh mode 
was added, also because of potential effect on the control system. Relative 
motion of fluid through the tank sump was introduced to facilitate later coup- 
ling with the feedlines. 

The structural idealization used in this study was developed during the 
Phase B Space Shuttle study by McDonnell Douglas Astronautics Company- East 
for booster Configuration 14 and high cross-range Orbiter 0050. The lumped 
mass and beam element idealization used is illustrated in fig. 2. Three 
degrees of freedom, axial, plunge, and pitch rotation were allowed at each 
mass. The fluid mass in each tank was elastically connected to the tank aft 
bulkhead so that its frequency corresponded to the first tank mode frequency 
predicted by a simplified hydroelastic model. The Phase B fully reusable 
baseline configuration, booster Model 256-20B, is illustrated in fig. 3. The 
existing data were extrapolated to the baseline design by scaling the Model 14 
length approximately 15 percent to make it conform with the Model 20B length. 
Natural frequencies and modal deflections remained unaltered. The natural 
frequencies for the model are presented in table 1. Modal damping of 
one percent was assumed for all modes, and all modes were mass normalized. 
Rigid-body axial, plunge, and pitch rotation modes were generated based on 
the mass properties generated during Phase B studies. 

The aerodynamic force at the center of pressure was determined from the 
relationship: 


REF 


q C l (M c 


) a 


( 1 ) 


where F ^ is the aerodynamic lift force, Srejt is the reference area, q is the 
dynamic pressure, C ^ is the lift coefficient gradient, M^ is the Mach number, 

and a is the pitch angle of attack. Numerical values for these parameters, and 
the location of the center of pressure were generated during Phase B studies. 
Aerodynamic damping forces associated with pitch rate ranged from 0. 03 to 
11.0 percent of critical damping. The fraction of critical damping was con- 
servatively set at 0. 01 percent, f, = 0. 0001. No other aerodynamic forces 
were considered relevant to the study. The aerodynamic forces of equation (1) 
are incorporated in the mathematical model as additions to the rigid-body 
partition (upper left) of the damping and stiffness matrices; 
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Figure 2. Space Shuttle Structural Idealization 
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Figure 3. Space Shuttle Phase B Baseline Configuration 




TABLE 1. NATURAL FREQUENCIES OF SPACE SHUTTLE STRUCTURE 
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Rigid-body axial, plunge, and pitch rotation structural modes (SI, S2, S3) omitted 
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and where r|g i» r ls2 ,r IS3 are the r igid-body axial, plunge, and pitch modal 
amplitudes;^ 4>© 3 are the plunge and pitch deflections of the rigid-body 
plunge and pitch mode shapes, and f is the magnitude of the distance from 
center of mass to the center of pressure. 
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An attempt was made to incorporate slosh mass because of potential 
control system coupling. Only the booster LOX slosh mass was incorporated 
since it is much larger than the LH 2 or orbiter slosh masses. It can be 
shown that the addition of a relative slosh degree of freedom, can 

accomplished by augmenting the structural equations as follows: 


° | V^SL ^ ZT 
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0 | 2 ^sl w sl 
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f ° 1 
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SL J 


q SL \ 


[ 0 1 



where Mg L is the slosh mass, £,sl is the slosh fraction of critical damping, 
W SL * s slosh frequency, Tg is the structure modal amplitudes, and <}> 
is the lateral structural mode shapes at the attach point. Slosh parameters 
were generated during Phase B studies. The slosh mass was attached at the 
bottom of the LOX tank since mode shapes at that point were already being 
used for feedline coupling. This choice of attach point led to a minor control 
instability. For reasons discussed in the baseline stability investigation 
section, slosh was omitted in subsequent runs. 

The structural modes as provided assumed no fluid motion out of the 

LOX tank; it was added since coupling the feedlines required that degree of 

freedom. The displacement of the tank mode degree of freedom, X^, is 

approximately the same as the motion of the fluid center of mass (c. m. ) with 

no flow. The total motion of the c. m. , is the sum of the tank mode 

motion and an additional motion of the fluid c. m. due to a small flow out of 

the tank, X For small relative flow, these relationships have the form: 

K 


- X L + X R 
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where 
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where is the relative mass flow, M. is the LOX tank fluid mass, lu is 
the distance from the c. m. to liquid surtace, 2 ApT is the area of the two 
feedlines, p is the density, is the LOX tank volume, qp-j, is the fluid 
displacement at the top of the feedline, X^ is the axial displacement at the 
tank bottom, Aj^ is an effective area ratio defined by the equation, and 4>XL’ 
4 > X t are the ax i a t structural mode shapes of the tank mode, and the tank bot- 
tom degrees of freedom. The tank fluid mass is removed from the structural 
coordinate and added to the total motion coordinate, and the spring due to 
ullage gas, Ky, is added. The structural equations are augmented as follows: 
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Expressing in terms of and X T yields: 
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where Qpp is the associated force, V q is the gas specific heat ratio, Pq the 
ullage gas pressure, Aq the tank area at the liquid surface, and Vq the 
ullage gas volume. Comparison of the stiffness additions, KyA^, with the 
accumulator and other system stiffnesses showed them to be two orders of 
magnitude smaller. Since the springs are effectively parallel, the ullage 
spring can be neglected. Transforming the mass additions to structural 
modal coordinates yields: 
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The equations of motion are augmented as indicated in equation (2) for 
aerodynamics, equation (3) for slosh, and equation (6) for relative flow. 

After incorporating these enhancements to the structural model, the equations 
of motion have the form: (double dashes indicate the compression of r|g in the 

fi r st te rm) 
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where the H's are the structural modal forces, I is an identity matrix, and 
’’lg, -^g, ^>g are the flexible modal amplitudes, frequencies, and dampings. 

Defining qg as the vector r)g, q the equation (7) may be summarized as: 
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F eedline 


Feedline dynamics are of primary importance to the pogo system. The 
feedline design (fig. 4) was idealized as illustrated in fig. 5. Because of 
potential interaction of the lateral runs with the structure, primary consider- 
ation was given to maintaining the effective lateral run mass and the total 
system mass. Only one of the two LOX lines was modeled due to symmetry. 

A detailed fluid dynamic model was developed of the feedline segments. 
The feedline model reflects the inertial and compliant characteristics of the 
fluid and the radial compliance of the line. The model does not include the iner- 
tial or compliant properties of the line itself, which are negligible compared 
with the structure to which it is attached. A compatibility relationship is used to 
couple the segments. The equations are transformed to relative modal coordi- 
nates. An asymmetric equation form is adopted because it facilitates coupling 
and provides pressures as explicit state variables. The form of the last coordi- 
nate is changed to engine weight flow to facilitate coupling to the engine. 

The feedline was modeled as 17 segments, each of which was modeled as 
a discrete mass and spring system. The fluid mass in each segment 
was divided equally among uniformly spaced discrete masses; the dis- 
crete spring rate, which includes the line circumferential elasticity was 
= (N - 1) E t p Aj-/ L (d p + E t) where N is the number of discrete masses, 

L is the segment length, p is the liquid bulk modulus, and E, t, d, Ap are 
the pipe's Young's modulus, thickness, diameter, and area. Modeling data 
for the feedline segments identified in fig. 5 are presented in table 2. 

V ' ' - - ... . . . 

Every segment has a mass at each end. To couple the feedline segments, 
a compatibility relationship is written at each joint stating the flow into the 
joint equals the flow out. The last mass of the upstream segment is then 
eliminated using the compatibility relationship. As an example of a joint 
with more than two segments, the relationship for C-D-P-Q is: 


A F, C (q F, C58 “ q R3 ^ " A F, D (q F, D1 " q R3* 

+ A F> p (q F> P1 - q R3 ) + a F; q (q F> - q R3 ) 
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TABLE 2. SPACE SHUTTLE LOX FEEDLINE MODELING DATA 
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or 


a f,d _ , a f,p 

q F, C58 A r „ q F,Dl + A 

I , X? , 


. F, Q 

q F , P 1 A f>c q F , Ql 


A F.D + A F.P + a f. q\ 

1 " XTT q R3 


F, C 


(9) 


where Ap ^ is the area of segment i, qp — is the displacement of the jth 
discrete mass in segment i, and q^- are the displacements of the feedline pipe 
due to structural motion at reaction coordinate i. There are five reaction 
coordinates: axial and lateral at the tank bottom, axial at the intersection of 
segments C, D, P, and Q, and axial and lateral at the lateral manifold, 

E-F-G. Deflections due to structural rotations were considered small in 
comparison to lineal deflections and were ignored. The feedline motions at 
the six engines were assumed to behave identically. This follows from sym- 
metry for pairs I-K, M-O, and P-Q, and is consistent with the small rotation 
assumption for the remaining groupings, I-K versus M-O, M-O versus P-Q, 
and I-K versus P-Q. That is, the torque due to thrust differential was 
neglected. 

To verify the assumption of no thrust differential torque, the feedline 
mode involving fluid motion between I-K and P-Q was isolated and coupled to 
the closest structural mode. The structural frequency was shifted to maxi- 
mize the coupling. Both analytical and numerical solutions revealed a highly 
stable system, verifying the assumption. If the accumulator compliance is 
eliminated, the feedline mode is raised beyond the frequency range of interest. 
With these assumptions, the bottom masses of segments I, K, M, O, P, 
and Q are coupled through six identical springs representing accumulator and 
cavitation compliance to a single coordinate, q^-, representing the average 
inertial displacement of fluid into the engines. 

The mass and stiffness matrix for the feedline segments were assembled 
as five reaction coordinates q^ through q^ 5 > the engine fluid displacement 
q^-, and 154 discrete coordinates, q^ through qj- Q ]4 • A similarity 
transformation with a matrix based on the eleven joint compatibility equations 
was used to couple the segments and reduce the system to qp^ through q^^, 
q and the 143 independent discrete feedline coordinates. The accumulator- 
cavitation springs were added to the stiffness matrix at this point. Terms in 
the first six rows and columns of the mass matrix were neglected. This is 
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equivalent to modifying the segment model such that the last mass in a seg- 
ment is zero, and augmenting the first mass in the next segment. After these 
operations, the feedline equations had the form: 



The normal modes of the cantilevered system, q Rl - ' • • q R5 " q C " °> 
were calculated. Equation (10) was transformed to relative modal coordinates, 
and the modal equations were premultiplied by the feedline mode shapes trans- 
posed, 4> ^ [see equation (11)]. 

Letting 

q F = q F " K FF K FRC q RC 


where 



yields 
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where Wp, ^p, Hp are the feedline natural frequencies, modal amplitudes, 
and generalized forces and where modal damping, £p = 0.01, was introduced 
in the equations. It may be noted that - Kj^£.p Kp?p Kpp^., the 

reduced reaction stiffness, is zero since the fluid is free to move in response 
to a reaction displacement. This unsymmetric formulation retains the reac- 
tion forces explicitly. 

• 

The weight flow through tins engine, W, is determined from the inertial 
displacement as Pg Ap (qp - Xp) where p is the density, g is the accelera- 
tion of gravity, Ap is the feedline area at the engine, and Xp is the engine 
axial displacement. With this substitution the coordinate transformation and 
the feedline equations of motion of equation (11) assume the form: 
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Control 


The Space Shuttle pitch plane thrust vector control system was used in 
this study, A block diagram of this control system is shown in fig. 6. The 
commanded and actual engine gimbal deflection are assumed identical since 
actuator resonances are above the frequency range of interest. However 
inertial effects of the gimbaled engine, sometimes referred to as "tail wags 
dog", were within the range of interest. The "tail wags dog" equations are 
developed and incorporated in a matrix format with the control equations from 
the block diagram. 

The mass matrix contributions of the engine were removed from the 
gimbal plane structural coordinates, Zq, and 9 q and assigned to the center 
of mass (c.m. ) engine coordinates, Zg and 9g. The additions to the struc- 
tural equations have the form: 
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where Mg and Ig are the engine mass and inertia, F C. is the control force, 

F a and M a are the actuator force and moment, and \ is the distance from 
c.m. to gimbal plane. The engine coordinates are related to the commanded 
gimbal angle, 6 , by: 0g = 0 q +6 , and Zg = Z ^ - \(0q + 6). Equation (13) is 
transformed to relative coordinates: ’ 
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where &ZG ant ^ Fqq are ne ^ lateral force and moment on the gimbal, 

T is the total thrust, 12 x 2.45 x 10° Newtons (12 x 550, 000 lb), and where it 
has been noted that internal actuator forces cancel and that the lateral con- 
trol force is approximated by the thrust times the angle (in radians) for 
small gimbal angles. The dynamic equation in 6 is dropped since the actual 
angle is assumed equal to the commanded angle. The quantities Mg\ and 
Ij- + w6 re computed indirectly from the tail wags dog frequencies, 

which are defined as the frequency at which the modal force is zero. 

If the gimbal feedback loop is reduced, the control system may be 
expressed as a matrix polynomial in the Laplace operator, s, and the dummy 
control variables, A, B, C, identified in fig. 6. 
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where OrQ’ ®ARU are rotation measured at the rate gyro and the 

attitude reference unit. 
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Engine 


Engine dynamics are of primary importance to the pogo system. Space 
Shuttle main engine dynamic characteristics were obtained from North 
American Rockwell Corporation, Rocketdyne Division. This engine model was 
of a preliminary nature and neglected cavitation effects, generally important 
in pump dynamic gains. Measurements on existing engines and other candi- 
date shuttle engine models show substantially greater gains. For these 
reasons the engine model was a primary candidate for parametric variation. 
As provided, the transfer function from LOX section pressure, Pg, to 
chamber pressure, Pq and LOX weight flow, W, are: 


and 



where 


M- 


6.4516 x 10 


A ** 

• 4 m 
sec 



The fuel side transfer functions provided by the engine manufacturer 
were at least two orders of magnitude smaller and were neglected. Thrust, 

T, is obtained from chamber pressure by multiplying by the effective area, 

0. 1181m (183.0 in. ). The transfer function P(g/Pg is presented as a func- 

tion of frequency in fig. 7 and fig. 8. It may be noted that a gain of at least 
0. 6 18 is required for the thrust oscillation to exceed the pressure times area 
oscillatory load acting down on the engine. Although the engine never reaches 
that level, other candidate engines exceed it. 

The axial force, F^q. * s sum the thrust and the product of the 
suction pressure, Pg, ancl the line area, 12 x Ajg, so the equation may be 
written 
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Figure 7. Space Shuttle Main Engine Gain Versus Frequency 




Figure 8. Space Shuttle Main Engine Phase Versus Frequency 
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where and K-gg, the coefficient matrices of the differential 

equation, denote the s^, s, and constant portions of the matrix polynomial 
equation. 


Coupled System 

The individual structural, feedline, engine, and control models 'are 
combined in one system of equations. To accomplish this, a global state 
vector was selected, and the feedline and structure are transformed to the 
state space. The feedline and structure are coupled and the control system 
is expressed in the global coordinates. Combining the engine, control, and 
coupled structure-feedline systems gives the final state equations of motion. 

The global coordinate system, or state space vector, was formed as 
follows: structural generalized coordinates, qg, feedline modes, t)p, J_,OX 

tank pressure, P^, engine suction pressure, Pg, engine weight flow, W, 
total thrust, T, gimbal rotation, 6 , control variable A, control variable B. 
This is a mixed state vector including both displacement and force type ele- 
ments. The formulation results in some dependent equations which were 
eliminated automatically by the computer algorithm. 


The feedline equations are first transformed to global coordinates. The 
engine and reaction coordinates are then expressed in terms of the appropri- 
ate rows of the composite structural mode shape matrix, 4>gp. which include 
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rigid, flexible, and slosh modes. The coordinate transformation and the 
equations of motion (see equation (12)) assume the form: 

if 

| q RE } = [^SRc] { q S ] 
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The force on the engine fluid displacement, Qr-, is the negative of the 
suction pressure, -Pg , times the area, 6 x Ag. The only source of modal 
force on the feedline, Hp, is the tank pressure, Pp times the area of the 
feedline at the tank, A- T , times the value of the feedline mode shape at the 
top of the line. Making these substitutions and transforming to modal struc- 
tural force yields: 


Let 



where 4 >sr is the first five rows of ^SRC* 
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Similarly, the structural equations are converted to global coordinates, 
using the first row of the coordinate transform in equation (18). 
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( 20 ) 


The generalized force on the tank fluid, Qpx> i s negative of the 
pressure, -Pt, times the area of the -two feedlines, 2 x App. The general- 
ized force on the structure has two sources, the negative of the reaction 
forces on the two feedlines, [equation (19)] and the gimbal forces. Combin- 
ing the feedline equations (19) with the structural equation (20) yields 
equation (21). 

Let 

Q S = ' 2<<> SR Q R + ^SXG F XG + *S0G F 0G + ^SZG F ZG 
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Then 




where «|>SXG> $S9G> $SZG are the axial, rotational, and lateral structural 
mode shapes at the gimbal. 
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Since the coefficients of and W are zero, the system can # be 
transformed from a state vector involving W to one involving W by moving 
the associated column of the mass matrix to the damping matrix. If this is 
done and the resulting matrices are denoted Mg, Dg, Kg, the equations may 
be written: 



( 22 ) 


The control system, equation (15), is expressed in global coordinates by 
noting that the rotation at the rate gyro and the attitude reference unit can be 
expressed in terms of the associated rows of the structural mode shape 
matrix, «}> SRG and 4» SARU> as 


e 

e' 


RG 

ARU 


*SRG _ 
ARU 




So 


0 

0 

0 


_sK e <,, SRG 


K e*SARU 


ojojo 

-i r- 

ojojo 


oio[ 0 
4 !— 


j 


i 


o; oi o 

_L l_l 

o|o!o 

I i I 


s 2 (I e +M e \ 2 ) 


-s 2 M e \-T 

016 s 2 - s -45. 
0 
0 


0 

0 


.24s + 15. 
.032s + 1. 
0 


0 

0 

-.24s - 15. 
0 

. 1 6 s + 1 . 


W 


< 

f xg 


6 G 


ZG 


6 

A 


0 

0 

0’ 

0 

la 


or 


K] ("} + [ D c] {*} + [ K c] f x l ■ {°} 


(23) 


where X denotes the state vector and.M^, D^, Kq the coefficients of the 
equivalent differential equation, denote the s2, s, and constant terms of the 
matrix polynomial equation. 

Combining these control equations (23) with the coupled feedline- structure 
equations (22) and the engine equations (17) yields: 
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The equations of motion for the two uncoupled systems, pogo alone and 
control alone, can be accommodated by dropping the rows and columns asso- 
ciated with $ and T, or with F Q „, F,_,_,, 6, A, and B, respectively. 


STABILITY INVESTIGATION- - 


The mathematical model developed in the previous section was used to 
investigate the nature of pogo-control interaction. The model was investi- 
gated in four configurations: control system only active, engine system only 

active, both systems active, and neither system active. This section dis- 
cusses the method selected for stability investigation and the selection of a 



baseline configuration. The sensitivity of the interaction and stability to par- 
ametric variations was investigated. The onLy deviation from a linear rela- 
tionship between the eigenvalues of the four configurations was observed in 
the region of very close roots. This interaction was further investigated 
through multiple parametric variations. The following symbols have been 
used to identify the curves in the figures: A, aero mode (the rigid-body pitch 

mode), Si, ith structural mode (the first three of which are rigid body), and 
Fi, the ith feedline mode. 


Method 


The system stability was evaluated on the basis of closed-loop roots. 
The type of linear relationship expected is discussed, and an interaction 
coefficient is defined to evaluate the deviation from linearity, which will be 
referred to as interaction. 


Stability characteristics were evaluated on the basis of the fraction of 

critical damping associated with the system eigenvalues. The eigenvalues 

st 


were determined by assuming an exponential solution, X - Xe° , to the 


The roots of the characteristic equation, 
provide the eigenvalues, s the non-trivial 


homogeneous system equations, 
determinant (s2 M + sD + K) =_0, 

solutions of (s? M + s^D + K)X^ = 0 provide the associated eigenvectors, X^. 

If a^ and ib^ are the real and complex parts of the ith eigenvalue, s^, the 
undamped natural frequency, oj and the fraction of critical damping, are 
defined by 


and 



(25) 


A negative fraction of critical damping indicates an unstable system; 
positive values indicate stability. 

The eigenvalues of the structural system are dominated by the imaginary 
part for a lightly damped structure. A small perturbation to the system 
equations will produce a small variation in the eigenvalues. The change in 
frequency is small. However, since the real part of the eigenvalue is already 
small, the change in damping may be relatively large. If two small perturba- 
tions are made to the system equations, the variation in the eigenvalues may 
be approximated by the sum of the two perturbations applied independently. 
This corresponds to retaining only the linear terms in a power series 
expansion of the characteristic equation. Since the frequency is nearly 
constant, the variation in damping is also approximated by the sum of 
the variations associated with applying the perturbations independently. 

In particular, if the fractions of critical damping in any particular mode 
are: £,q for the structural model alone, t,p for the structure with the 
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engine system, £, q for the structure with control system, and £pc for 
the combined system, the following relationship should be approximately 
true : 


^>PC ~ *= o + + 


or 


(26) 


ipc ~ ^ p + ^ O 


To provide a measure of the accuracy of this approximation, an inter- 
action coefficient, IC, is defined for each mode as: 


ic = ^P + t=c ~ ^PC 


IC 


f, p + 4 c ' & PC 


1 


(27) 


A positive interaction coefficient indicates a less stable combined system 
than would be expected by the linear combination hypothesis; a negative coef- 
ficient, a more stable combined system. For those modes in which the inter- 
action coefficient is small compared to one, the linear combination hypothesis 
is supported. 


Baseline 


To improve the computational efficiency, or to make the system more 
relevant to the interaction problem, a baseline configuration was selected 
about which subsequent parametric variations would be made. The slosh 
mode was eliminated, the accumulators were removed, the number of modes 
retained was reduced, and the 25 percent burn time was selected. 

The eigenvalues of the coupled system as developed in the analytical 
model had a marginally unstable slosh mode with respect to the control sys- 
tem, The source of the instability was traced to locating the slosh degree of 
freedom too far aft. To evaluate the importance of slosh, the system eigen- 
values were recomputed without slosh. A comparison of the eigenvalues with 
and without slosh showed no significant effect. The 25 percent burn time 
comparison is presented in fig. 9. The slosh mode was therefore deleted 
from subsequent runs. 
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Figure 9. Effect of Slo*h on System Roots 











The coupled system now had no instabilities. To evaluate the effect of 
accumulators, the system eigenvalues were computed with the accumulators 
removed. The system was still stable as illustrated in fig. 10. It therefore 
appears that with the preliminary engine model, no accumulators would be 
required. To enhance the pogo participation for this study, subsequent runs 
were made with no accumulators. 

The natural frequencies and those fractions of critical damping differing 
significantly from the input values of one percent are presented as a function 
of burn time in fig. 11 and 12, respectively. The lines trace a system eigen- 
value, and the test point symbols identify the predominant input mode in the 
system eigenvector. The effect of the control system is most pronounced at 
25 percent burn, which was therefore selected for the baseline. 

The preliminary version of the coupled system included 18 structural 
(to 15. 2 Hz) and 12 feedline (to 44 Hz) modes. To improve computational 
efficiency, this was reduced to 14 structural (to 9. 8 Hz) and two feedline (to 
7.4 Hz) modes. The resulting eigenvalues were compared with the more 
complete set. Both frequency and damping were found to be the same within 
an error of 2 percent. The reduced set of modes was therefore used in all 
subsequent calculations. 

The baseline system was 25 percent burn time with no slosh and no 
accumulators. The system eigenvalues are presented in fig. 13 for the 
coupled system and for the pogo and control systems individually. As 
expected, the coupled system deviation from the input one percent damping is 
closely approximated by the sum of the individual system deviations. A 
numerical calculation confirms this: all interaction coefficients were 0. 03 or 
less. Thus the linear approximation is accurate to within three percent 
(0. 03 percent error in fraction of critical damping) for the baseline system. 


Sensitivity 


It is recognized that the Space Shuttle configuration and the input par- 
ameters are preliminary in nature. Reasonable variations were made to the 
input parameters to investigate both the stability and the interaction. Par- 
ametric variations— were made to engine gain and damping, - control gain and 
frequency, feedline and structural damping and gain, and structural and feed- 
line frequencies. 

As previously noted, a pump inlet pressure perturbation does not pro- 
duce a net positive gimbal force for this engine, although it does for other 
candidate Space Shuttle engines. Typical gains (P^/Pg) for other engines are 
three to five times greater than the preliminary values for this engine. The 
coupled system damping is plotted as a function of engine gains from one to 
ten times nominal in fig. 14. The second feedline mode becomes unstable, 
i. e. , pogos, at a little over five times the nominal gain #> This indicates 
accumulators may be required if the engine performs more typically than pre- 
liminary models indicate. The coupled system damping for engine dampings 




Figure 10. Effect of Feedline Accumulators on System Roots 
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Figure 13. Coupled Versus Uncoupled System Roots 
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Figure 14. Coupled System Stability Versus Engine Gain 
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from 0. 1 to 10 times nominal are presented in fig. 15. The primary effect 
of engine damping is to increase the damping in the first feedline mode, thus 
reducing the likelihood of pogo. The interaction coefficient associated with 
varied engine models were all at least as small as for the baseline. 

The characteristics of the control system are primarily dictated by the 
gimbal angle feedback loop. Recalling fig. 6 the block has the form: 

3/(1 + 0. 016 s) = X /(I + e s) 


The control system gain is proportional to 1/X, and the rolloff frequency is 
proportional to 1/e . System damping is presented as a function of control 
system gain in fig. 16. Reducing the gain substantially reduces the stability 
of the aero mode, A, (rigid-body pitch) while increasing it reduces the sta- 
bility of the second flexible mode, S5. Fig. 17 presents system damping as 
a function of cutoff frequency. Reducing the cutoff frequency tends to destabi- 
lize several flexible modes, while increasing it has little effect. It therefore 
appears that the control system design is adequate. Modifying the control 
design affects the coupled system in the same way as the control system 
alone. Stability with respect to control design can therefore be evaluated 
based on the control system alone. Again, interaction coefficients were 
almost as small as those for the baseline, producing less than 6 percent 
error. 

The potential for substantial design modification made substantial varia- 
tions in structural and feedline parameters desirable. The feedline damping 
was varied from half to twice nominal (see fig. 18). The associated system 
roots varied correspondingly. Close structural roots were affected in a sim- 
ilar manner but to a lesser extent. Similarly, the structural damping was 
varied from half to five times nominal. The coupled roots varied correspond- 
ingly as illustrated in fig. 19. Reducing structural damping produces a con- 
trol instability with the second flexible mode, which is normally only 
0. 27 percent damped due to control coupling. The feedline roots were largely 
unaffected and therefore not plotted. Interaction coefficients were all at least 
as small as for the baseline. 

The system stability is affected not only by the input damping but by the 
modal mass of the input modes. The effective modal mass was varied by 
scaling the mode shapes, and thus altering the modal gain. The effect of 
varying feedline gain is illustrated in fig. 20. The predominant effect is to 
enhance the engine damping in the first feedline mode. Other modes are not 
strongly affected. The interactions remain as small as for the baseline. 

The effect of varying structural gain is illustrated in fig. 21. The effect 
is similar to varying control gain. The aero mode, A, system damping 
drops with the gain, and the second flexible structural mode, S5, destabilizes 
as the gain is increased. These effects are due to using a control system 
which is inappropriate for the modified structural system and would be elim- 
inated in the control design cycle. The first significant interaction, a 21 per- 
cent deviation from linearity, occurred at four times the nominal structural 
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Figure 15. Coupled System Stability Venus Engine Flow 
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Figure 17. Coupled System Stability Versus Control System Cutoff Frequency 
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Figure 19. Coupled System Frequencies and Dampings as a Function of Structural Damping 
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Figure 21. Coupled System Stability Versus Structural Gain 
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gain in the third flexible structure mode and the first feedline mode. This 
interaction is attributed to increased coupling between the two modes. The 
nature of the interaction is best interpreted after examining the sensitivity 
to shifts in input frequency. 

The root behavior when the structural frequencies are scaled is illustrated 
in fig. 22. A control instability will develop in the eighth flexible mode, Sll, 
as the structural frequencies are decreased. As the frequency drops, the 
predominantly first feedline root, FI, changes character and becomes pre- 
dominately the fourth flexible structural mode, S7, and vice versa. In the 
region of transition, an interaction (not illustrated) of 20 percent occurs in 
the two modes, S7 and FI. Interaction again occurs only for the strongly 
coupled modes. The effects of varying the feedline frequency are illustrated 
in fig. 23. To maximize the interaction, the feedline frequencies were set at 
structural frequencies; FI = S5, S6 and F2 = Sll, S12, S13. The character 
of several of the system roots changes as the input feedline frequencies sweep 
through the structural frequencies. As expected, substantial interactions 
occurred at several points where roots were strongly coupled. A 17 percent 
interaction was observed at FI = S5, and a 50 percent interaction was observed 
at FI = S6. 

Two types of instabilities were encountered during the sensitivity 
investigation. In one the coupled system instability results from a control 
instability, in the other from a pogo instability. Control instabilities from 
varying control parameters or structural parameters all occurred when the 
control system was inappropriate for the associated structure. Pogo 
instabilities resulting from gains more representative of existing engines 
would be eliminated by accumulator design. Frequency shifts resulting in 
strongly coupled, closely spaced roots yielded strong interactions; in all 
other cases the linear combination of individual dampings was found to be 
very accurate. To investigate the region of interaction in more detail, the 
multi-parameter runs described in the next section were made. 


Interaction 


The nature of the deviation from linearity, which is referred to as 
interaction, in the region of significant coupling between closely spaced roots 
was investigated. The root behavior was investigated in detail as the feedline 
frequency was slowly swept through the structural frequency. The engine 
model was simultaneously set at four times the gain and one-half of the damp- 
ing of the nominal engine. This enhances the pogo contribution and makes it 
more representative of typical engines. 

Two regions were investigated in detail, the fir st and second feedline mode s 
in the vicinity of the second and eighth flexible structural modes respectively. 
The individual and coupled system root behavior in the vicinity of FI = S5 
and F2 = Sll are illustrated in fig. 24a and fig. 24b, respectively. Fig. 24a 
illustrates the fact that a stable control system and a stable pogo system can 
be combined to produce an unstable coupled system. This instability would be 
predicted on the basis of a linear combination outside the transition region. 

It may be noted from the figures that the damping crossing occurs at different 
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Figure 22. Coupled System Roots as a Function of Structural Frequency 
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Figure 24B. Coupled and Uncoupled System Stability Versus Feedline Frequency 
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points for each individual system and for the coupled system. The mode shape 
for each of the root pairs is undergoing substantial changes in response to 
small variations in the feedline frequency, moreover the transition is occurr- 
ing at different input values for each of the individual systems and the coupled 
system. It is not surprising that the linearity assumption breaks down in that 
case . 

The interaction coefficients associated with fig. 24a are presented in 
fig. 25. The interactions for the root pair are equal in magnitude but oppo- 
site in sign and appear to be correlated with the frequency crossing in the 
control system, however further study would be required to fully understand 
the behavior of the interaction coefficients. The linear combination of the 
envelopes of the separate systems gives the envelope of the coupled system 
roots. The envelope of the root behavior is meant to imply the smooth curve 
extrapolated from the region on either side of the associated root crossing. 
Fig. 26 illustrates the extrapolated control and pogo system curves and the 
linear prediction for the coupled system. The actual coupled root behavior 
is depicted by the broken line. 


CONCLUSIONS 


The nature of the pogo- control interaction was investigated using a 
coupled structural, feedline, control, and engine system model developed 
for a representative Space Shuttle. To evaluate the sensitivity of the inter- 
action and stability to design modification, the model was subjected to sub- 
stantial parametric variations. The stability of each system was determined 
by closed-loop eigenvalue analysis for four conditions: engine system only 
active, control system only active, both engine and control systems active, 
and neither engine nor control system active. Based on these investigations, 
the following conclusions may be made about vehicles with substantial asym- 
metries or lateral feedline runs: 

(1) The coupled pogo-control system may be unstable even though the 
pogo and control systems are separately stable. Major instabilities 
may exist even though conventional analysis of both the pogo and 
control system has shown them to be stable. 

(2) The coupled pogo-control system stability can be evaluated on the 

basis of the separate stabilities of the pogo and control systems as 
determined by conventional analysis techniques. Specifically, the 
coupled pogo-control system fraction of critical damping, for 

any mode is the sum of the damping with the pogo system only, £,p, 
plus the damping with the control system only, 1, q , less the structural 
damping with neither pogo nor control systems, £,q. 
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Figure 25. Interaction Coefficient Versus Feedline Frequency 
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Figure 26. Extrapolated System Stability Versus Feedline Frequency 
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(3) Although this relationship is imprecise for very close roots, it may 
be used to predict the upper and lower bounds for the root pair. The 
individual pogo and control system roots are extrapolated through 
the region of strong coupling from points outside the region. The 
upper and lower bounds for the coupled root pair are determined by 
applying the linear relationship to the extrapolated curves. 

On the basis of these conclusions it is evident that pogo- control coupling 
must be evaluated for vehicles with structural asymmetries or long lateral 
feedline runs. It is suggested that conventional analysis be performed on the 
separate pogo and control systems, and that the derived relationship be used 
to verify the stability of the coupled system. Closed-loop eigenvalue analysis 
was found to be a very effective technique, and is recommended for future 
stability investigations. Although the understanding of interaction is incom- 
plete, enough insight has been gained to be used effectively in the design pro- 
cess. It is therefore recommended that future effort be concentrated in the 
less well understood aspects of pogo such as asymmetric tank dynamics and 
engine dynamics. 


McDonnell Douglas Astronautics Company, 
Huntington Beach, California, 

May 9, 1972 
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